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Abstract 



This paper extends earlier work on the definition of Wannier functions for Bloch 
electrons in a magnetic field. Extensions to irrational as well as rational magnetic fields 
are defined, and their properties investigated. The results are used to give a generalisation 
of the Peierls effective Hamiltonian which is valid when the magnetic fiux per unit cell is 
close to any rational number. 
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1. Introduction 



Wannier functions are localised basis states which span a band of Bloch eigenfunc- 
tions [1]. The use of localised basis functions can be convenient both technically and 
conceptually, particularly when considering perturbations which are themselves spatially 
localised. There are difficulties in defining satisfactory Wannier states when a magnetic 
field is applied to the lattice. Firstly, the eigenfunctions are typically not Bloch states: in 
two dimensions the eigenfunctions are only Bloch states if the ratio (3 of the flux quantum 
to the magnetic flux per unit cell is a rational number [2,3] (in these cases I will write 
13 = p/q, where p and q are integers with no common divisor). Secondly, even when the 
magnetic field is rational in this sense, conventional Wannier states only have satisfactory 
localisation properties if a topological invariant (the Chern index) chararterising a Bloch 
band is equal to zero [4,5]. In a previous paper (reference [6]) it was shown how this latter 
difficulty could be overcome, for two dimensional lattices, in the case where the magentic 
fiux per unit cell is rational. In [6] I showed how to obtain a complete set of states which 
span a Bloch band, and which retain all of the useful properties of conventional Wannier 
functions. Two different definitions were examined, termed type I and type II Wannier 
functions. The definition of these states contains the Chern index M of the band, and in 
both cases they reduce to the conventional Wannier function when M — Q. 

The purpose of the present paper is twofold. The first objective is to show how the 
definition of Wannier functions can be usefully extended to irrational fields, despite the fact 
that Bloch bands do not exist in this case. Some of the results for type II Wannier functions 
are anticipated in earlier papers by the same author [7] , [8] (the latter in collaboration with 
R. J. Kay). These earlier papers discussed 'irrational' generalised Wannier functions for 
the special case of the 'phase space lattice Hamiltonian', a one dimensional model which 
reprenents many of the features of Bloch electrons in a magnetic field. The form of the 
Wannier functions of the full Hamiltonian for irrational magnetic fields is related to the 
results for the phase space lattice Hamiltonian, but the generalisations are not obvious. 
The derivation given here is also more satisfying in that it uses only minimal algebraic 
properties, and that results are obtained for both types of Wannier functions introduced 
in [6]. 

The second objective is to use the generalised Wannier functions to obtain a very 
general form of the Peierls effective Hamiltonian [9,10], in a form suitable for system- 
atic analysis. For simplicity we consider only a two-dimensional case where the electron 
is confined to a plane, and perturbed by a magnetic field in the perpendicular direction 
(with cartesian coordinate z). A comprehensive treatment of the three dimensional case 
introduces the complication of an additional commensurability parameter, but is straight- 
forward when the field is aligned with one of the crystal axes. The Peierls Hamiltonian is 
a one-dimensional effective Hamiltonian which describes the effect of a uniform magnetic 
field perturbing a band of Bloch states. If the dispersion relation is S{kx,ky), the Peierls 
effective Hamiltonian takes the form 

where Kj^ and Ky are generators of the magnetic translation operators, T(R) (These are 
defined in section 3; they were introduced in [11,12], and are discussed concisely in [6]). 
These satisfy 

where Aj are the basis vectors for the lattice. Many derivations of this relationship exist 
where the dispersion relation is that of the B = problem. This paper considers the case 



where the dispersion relation is that of the system with any rational magnetic field p/q, 
showing that the Peierls effective Hamiltonian is applicable in this case. It is shown that 
the commutator (1.2) is replaced by one which depends upon the Chern index: the general 
form of (1.2) is 

where 7 is another dimensionsionless parameter characterising the magnetic field. The 
value of 7 depends upon the value of the Chern integer M, and upon another integer N 
which satisfies 

qM + pN = l. (1.4) 
The dimensionless effective magnetic field 7 is 

^ M + NP ' ^ ' ^ 

The expression (1.5) can be surmised from results obtained previously for the phase space 
lattice Hamiltonian [7,8]. The derivation presented here indicates how the effective Hamil- 
tonian can be obtained for the full Hamiltonian, rather than a one-dimensional model. 
This issue has also been considered by Chang and Niu [13], who also discussed an heuris- 
tic approach to determining the first order correction to the effective Hamiltonian. The 
method described here allows a systematic development of the effective Hamiltonian, using 
similar techniques to those applied to the phase space lattice Hamiltonian in reference [14] . 
It also has the advantage that some of the complicated intermediate steps in the algebra 
of references [7] and [14] are given a more transparent interpretation. 

Sections 2 and 3 respectively summarise the essential definitions and principal results 
from [6] , and a representation of the Hamiltonian as a sum of magnetic translation opera- 
tors. The latter will be essential to the derivation of the general effective Hamiltonian. 

Section 4 describes the extension of the Wannier functions obtained in [6] to irrational 
magnetic fields, and a corresponding extension of the definition of Bloch states. The next 
four sections consider various properties of the generalised Wannier functions and Bloch 
states. Section 5 discusses the effect upon the Wannier functions of a transformation of 
the Bloch states. A 'gauge transformation' of the form 

|5(k))^|S'(k))=exp[i^(k)]|S(k)) (1.6) 

is applied to the Bloch states, with ^(k) a periodic function. The Wannier states derived 
from the gauge transformed Bloch states can be obtained from the original Bloch states by 
the action of an operator, which is obtained in section 5. Similarly, section 6 determines an 
operator acting on the Wannier states which is the image of a translation operator acting 
on the Bloch states. 

Section 7 computes the Dirac bracket of two generalised Bloch states, (i?'(k') |i?(k)), 
which will be required for determining matrix elements of the Hamiltonian. Section 8 
introduces some notational devices which simplify and illuminate the rather complex ex- 
pressions obtained earlier, showing how they can written in terms of translation operators 
with algebra analogous to that of the magnetic translation group. Finally, in section 9 
these results are used to obtain the general form for the Peierls efffective Hamiltonian. 



2. Summary of earlier results 



The purpose of this section is to present, for the convenience of the reader, a summary 
of some of the principal definitions and equations from the earher paper, reference [6]. 
The lattice vectors are written R = niAi + n2A2, and the reprocal lattice vectors are 
K = niai + ?i2a2, with a^.Aj = 2n6ij. 

The magnetic translation operators T'(R) introduced by Zak [11] and Brown [12] are 
of fundamental importance. They are a representation of the symmetry of the system: if 
R is a lattice vector, the T'(R) commutes with the Hamiltonian. The magnetic translation 
operators do not commute among themselves, and their composition rule can be written 
in the form 

>i (Ri A R2) 



T(Ri)T(R2) = exp 



P (Ai A A2) 



T(Ri+R2) (2.1) 



where P is the flux quantum divided by the magnetic flux per unit cell. The magnetic 
translation operators are discussed concisely in [6] . 

When conventional Wannier functions are defined, it is assumed that the Bloch states 
are periodic functions of the Bloch wavevector k, as well as being eigenfunctions of the 
lattice translation operators T{Ai), with eigenvalues exp[ik.Ai]. In the case where a 
rational magnetic field (with q/p flux quanta per unit cell) is applied, in general both of 
these conditions need to be modifled. The Bloch states are p fold degenerate, and their 
phase increases by 27rM on traversing the boundary of the unit cell. Throughout this 
paper, the following choice for the eigenvalue and periodicity conditions is preferred: 

f(Ai)|S(k)) = exp[ik.Ai]|S(k-ga2/p)) (2.2a) 

f (A2) |S(k)) = exp[ik. A2] |S(k)) (2.26) 
\B(k + ai/p)) = exp[iMk.A2] |S(k)) (2.2c) 

|S(k + a2)) = |B(k)) . (2.2d) 

Bloch states with their phases chosen to satisfy (2.2a) and (2.2b), and with degenerate 
states resolved so that (2.2c) is satisfied will be termed canonical Bloch states. Except 
when p = 1 and M = 0, these conditions depend upon the choice of lattice basis vectors 

The method for constructing the Wannier functions is based upon the following ob- 
servation: if the Bloch states are canonical, the state |C(k)) = T(— pMk.A2/27r)|-B(k)) 
is periodic on the Brillouin zone of the superlattice spanned by pAi, A2, and Wannier 
functions |x(I^)) obtained by integrating the state IC'(k)) with weight exp[ik.R]. In 
the case of standard Wannier functions, all of the Wannier states are obtained be applying 
translation operators to a single fundamental Wannier state. In the magnetic case, the full 
set of Wannier states is obtained by applying lattice translations to |A^| fundamental type 
I Wannier states, \x^) — |x(/^Ai)), /x = 0, .., |A^| — 1, where N satisfies (1.4). The relation 
between the Bloch and type I Wannier states is 

|5(k)) = exp[-ik.R]T(n2A2)T(niAi)f (^(k.Ai)A2) 



R=ni Ai-|-n2A2 

\N\ 



J2 exp[-ip//(k.Ai)]|x/.) . (2.3) 



A somewhat more natural representation of the Bloch states uses an alternative set of 
fundamental Wannier states: the type II Wannier states are defined by 

|JV|-1 

1*^/^) =nY1 exp[-27ri^/x7iV]T(-/.'Ai/7V)|x^.) . (2.4) 

One advantage of using the type II Wannier states is that upon expanding the Bloch states 
in terms of the states, the summation over /j, no longer depends upon k: the Bloch 
states are given in terms of the type II states by the relation 

|7V|-1 

\B{k)) = exp[-ik.R] ^ exp[27rini///A^] 

R=niAi/iV+n2A2 M=0 

xf(n2A2)f(niAi/iV)f (^(k.Ai)A2)|</)p) . (2.5) 

The other advantage of the type II Wannier states is that their transformations under a 
change of lattice basis vectors are simpler [6] . 



3. The Hamiltonian in terms of translation operators 

Here the objective is to represent the Hamiltonian as a sum of magnetic translation op- 
erators: this will facilitate the construction of the effective Hamiltonian. The Hamiltonian 
is 

F(r) = y(r + R), R = niAi+?i2A2 (3.1) 

with the magnetic field generated by a linear vector potential, constructed using a matrix 
B with elements Biji 

A(r) = Br, V A A = Beg, B21 - B12 = B . (3.2) 

The magnetic translations T'(R) have a generator P = Piei + P2e2: 

T(R) = exp[-iP.R/n] (3.3) 

P = p - eB^f . (3.4) 
It will also be useful to define a set of conjugate generators P*: 

P* = p - eBf . (3.5) 

The generators Pi, P^ satisfy the commutation relations, where Sij is the antisymmetric 
symbol, with elements £11 = £22 = 0, £12 = —£21 = 1: 

[Pi,Pj] = -iehBsij (3.6) 
[P:,P*]=iehBeij (3.7) 



[Pi,Pj*]=0. (3.8) 

The coordinate vector f can be expressed in terms of the generators P and P*: from (3.4) 
and (3.5) it follows that P* — Pi = eBsijrj (where, from here until the end of section 3, 
repeated indices are summed over). This can be inverted to give 

ri = ^eij{Pj-P;). (3.9) 
The Hamiltonian can now be written 

H = ^P:p: + E exp[ik.r] (3.10) 

k 

where the k = niai + n2a2 are vectors in the reciprocal lattice, with basis vectors satisfying 
a^.Aj — 27T5ij. Expressing the r using (3.9), and using the fact that Pi and Pj commute, 
(3.10) can be written in the form 



H = ^P:p: + J2'^^ ^M^kiSijPj/eB] exp[-ik,eijP;/eB] 

k 

= T^PiPi + E f*{-nk*/eB) f{nk*/eB) (3.11) 

k 

where k* = k*ei and T*(R) are defined by 

k* = Eijkj, f *(R) = exp[-iP*.R/n] . (3.12) 

The Hamiltonian is therefore expressed in terms of a sum of magnetic translation operators, 
with operator- valued coefficients Vk^ 

H= J2 14 f{Kk*/eB) . (3.13) 

k=niai+n2a2 

The operators Vk commute with the magnetic translation operators, and are given by 

Vk = -3-5k,oP*.P* + V^k f*{-Kk:*/eB) . (3.14) 
Am 

It is desireable to express the vectors k* in terms of the real-space lattice basis vectors Ai, 
A2. The vectors corresponding to reciprocal lattice vectors are denoted by k*. Writing 
Aj = ^ij^j and a^ = aijBj, the matrices A — {Aij} and d — {a^} satisfy Ad^ = 2nl. It 
follows that ^ 

'^i = TT7TV^~^2iei - ^2262) 
det(^) 

•^2 = —--riAiiei + ^1262) . (3.15) 
det(^) 

Noting that the det(JL) is equal to the area A of the unit cell, the Hamiltonian (3.13) can 
then written 

00 00 / u \ 



where k = njaj. The elementary transformations associated with this representation of 
the Hamiltonian, Ti = T{—hA2/eBA) and T2 = T{hAi/eBA) therefore span a lattice 
which is aligned with the crystal lattice, but scaled by a dimensionless factor 

^ h flux quantum 

eBA flux per unit cell 



4. Extension to irrational magnetic fields 

When the number of flux quanta per unit cell is rational, the spectrum consists of Bloch 
bands for which Wannier functions have been defined. When the number of fiux quanta 
per unit cell is irrational, there are no Bloch bands and the spectrum is a Cantor set. It 
is however still possible to define useful sets of generalised Bloch states and corresponding 
Wannier functions. 

The expression giving the Bloch states in terms of the type II Wannier states will be 
generalised, by writing 

|JV|-1 

\B{k)) = exp[-ik.R] exp[27rini///iV] 

R=niAi/Ar+n2A2 M=0 

xT(n2A2) T(niAi/iV) T(M(k.Ai)A2/«) |</.^) (4.1) 

Straightforward application of the composition law (2.1) for magnetic translations to the 
form (4.1) shows that the generalised Bloch states satisfy a periodicity condition 

|S(k + Kai/27r)) = exp[iM(k.A2)] |5(k)) (4.2) 

provided that exp[27riMn//3A'"] exp[— i^n/A^] = 1 for all integer n. The latter condi- 
tion is used to determine allowed values for the constant k: this quantity must satisfy 
k(3 = 2'k{M + I3NJ) with J an integer. Equation (4.2) is a natural generalisation of the 
periodicity condition (2.2c). It is desireable to define the generalised Bloch states so that 
as (3 ^ p/q they converge to the Bloch eigenstates of the rational case with (3 — p/q. 
Setting J = 1 (and using (1.4)), n aproaches 27r/p as ^ p/q, which is consistent with 
(2.2c). The appropriate choice of the constant k, defining the dimension of the Brillouin 
zone is therefore 

K/3 = 27r(M + /3N) . (4.3) 

Systematic application of (2.1) shows that the states (4.1) also satisfy other conditions 
analogous to the standard Bloch states: collecting together the periodicity properties and 
the equations defining the effect of lattice vector translations, the generalised Bloch states 
satisfy the relations 

|5(k + Kai/27r)) = exp[iM(k.A2)] |5(k)) (4.4a) 

|S(k+a2)) = |S(k)) (4.46) 

T(Ai)|S(k)) =exp[i(k.Ai)]|S(k-a2//3)) (4.4c) 

T(A2)|i?(k)) = exp[i(k.A2)]|S(k)) . (4.4d) 

Equation (4.1) defined the generalised Bloch states in terms of type II Wannier states. 
Using the relation between the type I and type II Wannier functions given by (2.4), the 



(3.17) 



corresponding relation giving the generalised Bloch states in terms of type I Wannier 
functions is 

\B{k))= Yl exp[-ik.R]T(n2A2)T(niAi) 

R=niAi+n2 A2 

|Ar|-i 

xT(M (k.Ai)A2/Av) J2 exp[-27riM(k.Ai)/«]|xM) (4-5) 

On systematic application of (2.1) and (4.3), it is found that the states (4.5) satisfy the 
canonical Bloch state relations in the form (4.4a-d), for any states The relation 

between the type I and type II Wannier functions therefore remains valid in the irrational 
case. 

The generalised Bloch states lie in a Brillouin zone spanned by the reciprocal lattice 
vectors Kai/27r and sl2, with area = K\aii AsL2\/2n. Applying Born- von Karman bound- 
ary conditions, the density of states per unit area associated with the set of generalised 
Bloch states is ^k/47r^. The area of the real space unit cell, A = |Ai A A2I is equal to 
47r^/|ai A a2|. The density of generalised Bloch states states per unit area is therefore 



It will now be shown that this density of states is precisely what is required for them to 
form a complete set of states for a region of the spectrum bounded by two gaps. Stfeda [15] 
showed that the density of bulk states per unit area for region of the spectrum bounded 
by two gaps satisfies is related to the Hall coefficient a^y'- 

= e— . (4.7) 



The Hall coefficient is quantised in units of e^/h, and the Chern number M is the quantum 
number [5]: 

a^y = M—. (4.8) 

The density of states is clearly correct in the rational case. Using (4.3) with the relation 
(3 — h/eBA to differentiate (4.6) with respect to S, equation (4.7) reproduces (4.8). 
This shows that the variation of the density of generalised Bloch states with respect to 
magnetic field is precisely the same as that of the eigenstates. The generalised Bloch states 
are therefore a complete set provided that they are not linearly related. 



5. Images of gauge transformations 

The gauge transformations considered are of the form (1.6), in which the Bloch states 
are multiplied by a factor exp[i^(k)]. The cases of rational and irrational fields will be 
considered separately. 

5.1 Rational case 

In the rational case 6 satisfies 



^(k + ai/p) - 27rLi = ^(k) = ^(k + a2/p) - 27rL2 (5.1) 



with Li and L2 integers, so that the gauge transformation leaves the Bloch states in 
canonical form. Wannier functions may be defined for the gauge transformed states. These 
Wannier functions will be different from the original ones, and it is interesting to determine 
how the transformed Wannier functions may be obtained from the original ones directly. 
The calculation will be presented for the special case where 

exp[i6'(k)] = exp[ik.R*] (5.2) 

where 

R* = p(Li Ai + L2A2) (5.3) 

is a superlattice vector. More general transformations of the form ^(k) = k.R* + e6'(k), 
with 6'(k) periodic in ki and /c2, with period 27r/p, can also be treated for e -C 1 by Fourier 
expanding 6'(k). The type I Wannier functions of the gauge transformed Bloch states are 

IX'(R)) = £2 exp[ik.(R + R*)]T(^(k.Ai)A2)|i?(k)) 

= |x(R + R*)). (5.4) 

The fundamental type I Wannier functions, |x^), are a subset of the full set of Wannier 
states |x(R)), defined by \xij) — |x(p/^Ai)): previously the index jjL was restricted to the 
range /i G {0,..,|A^| — 1} but it is convenient to extend the definition by allowing jj. to 
take any integer value. The states |x(R')) ^i'^ obtained from the fundamental Wannier 
functions by the relation [6] 

\x{p{Nni + ;u)Ai + nsAs)) = f (risAs) f (mAi) |x^) . (5.5) 

It follows that the extended set of fundamental Wannier states satisfies 

Ix^+at) = T(Ai) Ix^) . (5.6) 

The transformation of the fundamental type I Wannier functions is therefore 

Ixp = T(pi:2A2)|xM+Li) • (5.7) 

The corresponding transformation of the type II Wannier states is obtained using (2.4) 
and its inverse relation as follows 



|iV|-l 

y = ]^ E exp[-27ri/x/x7iV]T(-/x'Ai/iV)|x;,,) 



/Lt'=0 

lATI-l liVI-1 

N 



= ^ E E exp[27ri(A-/x)/x7iV]exp[27riALi/iV] 

p'=0 A=0 

xf (VAi/iV)r(pL2A2)f ((//' + Li)Ai/iV)|(/)A) . (5.8) 

After combining the translation operators, the summations can be performed: only the 
term A = + qL2 contributes, giving the result 



10^) = exp 



27ri(/i + L2q')Li 
N 



f(pL2A2)f(LiAi/iV)|0^+,i,) 



f{LiAi/N + pL2A2)\4>^+qL,) . (5.9) 

This expression will be re-cast into a more transparent form in section 8. 
5.2 Irrational case 

Now consider the case of gauge transformations of the generalised Bloch states defined 
for irrational fields. In order to define a transformation of the Wannier functions, the gauge 
transformation must leave the Bloch states in canonical form. If /5 is irrational, equations 
(4.4b) and (4.4c) imply that a suitable gauge transformation cannot depend upon /c2- 
Linear gauge transformations analogous to (5.2) are therefore restricted to being of the 
form 

\B'{k)) = exp[27rifciLi/«;]|S(k)) . (5.10) 
Using (4.1), a Bloch state may be written in terms of the Wannier functions \(f)'^) as follows 

oo oo \N\-1 

\B'{k)) = ^ XI XI exp[-iA;i(ni - Li)/N] exp[-iA;2n2] exp[27ri//(ni - Li)/N] 

ni=— oon2 = — oo /n=0 

xT(n2A2)T((ni - L^)AjN)f{MkrA2/K) . (5.11) 
If the Wannier functions generating this state are 

10;) = exp[27ri/xLi/Ar]T(LiAi/Ar) (5.12) 

then (using (4.3)) it can be seen that \B'{k)) is related to the original Bloch state by (5.10). 
This result reduces to a special case of (5.9) in the case where ft = p/q. 

6. Images of translation operators acting upon Wannier states 

This section discusses the states 

f (r)|S(k)) , r = (3{v^A^ + v^A^) (6.1) 

with vi, V2 taking integer values. It will be demonstrated that they are generalised Bloch 
states of the form (4.1), generated by a set of Wannier functions \4>'^), fJ' = 0, ••, |-^| — 1- 
The transformation giving these Wannier states in terms the states |(/)^) which generate 
the original Bloch state |-B(k)] will be determined. This transformation may be regarded 
as the image of the operator T(r) acting on the Wannier functions. 

The wavevector k = (/ci, /C2) of the state (6.1) is shifted to {ki + A/ci, /C2), with Aki 
to be determined. Commuting the operator T'(r) to the right using (2.1) gives 

00 00 

T(r)|S(k))= ^ exp[-i(/ci + A/ci)ni/iV]exp[-i/c2n2]exp[27ri//ni/iV] 

ni = — 00 n2 = — 00 p=0 

xT(n2A2)T(niAi/Ar)T(M(A;i + AA;i)A2/«) 

exp[iA/cini/A^] exp[— 27riz/2?^i/-^] exp [27ri(A;i + ^AkiM)iyi/K] 

xfiPu^A^ + {(3v2 - MAki/n)A2)\(l>^) . (6.2) 



exp 



27ri(// + ^qL2)Li 



N 



This state is a generalised Bloch state if the product of the final two phase factors containing 
rii is unity. This occurs if Aki = 27rz/2- In this case the argument of the last translation 
operator simplifies, the multiplier of A2 becoming 27tNu2/3/ n. The state (6.1) is then in 
the form 

f(r)|i?(k))=exp[i0(k')]|S'(k')) (6.3) 
where |i?'(k)) is a generalised Bloch state with Wannier functions 10^, k' = k + ^'2ai, and 

«(k) = . (6.4) 

K 

The Wannier functions generating \B'(k)) are 

= exp[-27r2iMzyii/2/«]T(/?i/iAi + 27r7V/?zy2A2/«)|0^) . (6.5) 

The phase factor in (6.4) represents a gauge transformation of the type (5.10). Using 
(5.12), we may therfore write 

T(r)|S(k)) = |S"(k + iy2ai)) (6.6) 
where the Wannier states generating |-B(k)) are 

\(f)'P = eM^mMfiu^/N]f(Mu^A^/N)\cf)'^) . (6.7) 

The Wannier functions generating the Bloch states \B" (k + 1/23.1)) = T(r)\B(k)) can now 
be expressed in terms of the original Wannier states: 



{(Pa) = exp 



2'KiMfiui 
N 



,/ K/3 ^ 2nN/3 ^ \ ,^ , 



7. Dirac brackets of generalised Bloch states 

The objective is to evaluate the matrix element 

/(k,k') = (i?'(k')|i?(k)) 



(7.1) 



where the |-B(k)) and |5'(k)) are different generalised Bloch states for irrational magnetic 
fields. These Bloch states are generated by different type II Wannier states and 
\4>'^) respctively, using the expansion (4.1). The resulting expression will later be used to 

calculate matrix elements of the form (S(k') |T(r) |i?(k)), (where r = uiSli + 1^232), and 
hence matrix elements of the Hamiltonian, by writing |i?'(k + iy2^i)) = T{r)\B(k.)) . 
Using (4.1) and (2.1), and writing ki = k.Aj, the Dirac bracket is 



00 00 



00 00 |iVI-l|JV|-l 

'= j: e e e e e-p 

ni=-oo n'i=-oo n2 = -oo n^ = -oo M=0 m'=0 



i{k'2 - k2) 



X exp 



27r / 
i{{k[ - ki) + 27r(/x - /x') - — (n2 - na)) ( 



ni + n[ 
2N 



exp 



n2 + n'2 



ki + k[ 
2k 



(n'l - ni) 



X exp 



27ri f n + n' 



N \ 2 

, 1^/ ni -n[ 



(ni -n[) 



exp 



(n2 - 



M 



It is convenient to make changes of variable 

j = ni- n[ J = 
I — n2 — n'2 L = 



(7.2) 



ni + n'l 
2 

77-2 + 77-2 



(7.3) 



The summations in (7.2) will then run over integer values of L for even I, and over integer 
plus one half values of L for odd I, similarly for J and j. These sums are most conveneiently 
evaluated by decomposing them into four summations: 

j I J L 



00 00 00 



= E E E E 

n——oo n'=— 00 m=— 00 m'= — oo 



F{2n, 2m, n', m') + F{2n + l,m,n' + ^, m!) 



+F(2n, 2m + 1, n', m' + i) + F(2n + 1, 2m + 1, n' + |, m' + i) 
The function F(j, l,J,L) is of the form 

-P'(j, ^, J, L) = exp[iai(Z) J] exp[ia2-^-] f{j, 

where 



iV 

CK2 = ^2 ~ ^2 



(7.4) 



(7.5) 



(7.6) 



and 



7V|-1 |Af|-l 
H=0 ij,'=o 



2k 



2Ar 



exp 



k2 + k'2 



3 . M 



x(<^;,|T(^Ai + {l+-{k^- k[))A2)\cl>, 



K 



(7.7) 



Using (7.5), it is seen that the sums over J and L are easily evaluated using the Poisson 
summation formula in the form 



exp(iQ;n) = 27r 5 (a — 27rm) 



(7.8) 



m= — oo 



Using this formula, 

. 2 °o °o 00 00 
^(k',k) = ^ E E E E (-l)^"^^^+^^^^^<^(/^2 - - 27riV2) 

iVi=— 00 iV2=— 00 ni = — 00 n2 = — 00 



\N\-1 \N\-1 



2tt 



X ^ ^ S{ki -k[- 2tt{ii - ijl') + —77,2 - 2-kNNi) exp 

/x=0 /x'=0 ^ 



n2 



X exp 
Writing 



-27ri 



ki + k[ 11 + 11' 
2k 2 



ni 



+ (n^ + ^(/ci - fcl))A2)|0^) . (7.9) 



A/c = 27r(g-^) 



(7.10) 



and recalling (4.3), the values of ki — k[ for which the matrix element is non-zero may be 
written in two alternative forms: 



1 



ki - k[ = hAk + I2K = 27r(Li + -L2) 



(7.11) 



where Zi, I2 and Li, L2 are all integers. The argument of the second delta function in (7.9) 
can therefore be written in terms of AA; and k. Noting that 



d{hM) 



q N 
-p M 



= 1 



(7.12) 



it is seen that the sums over A^i, fi' and ?i2 in (7.9) may be replaced by a sum over the 
indices /i, I2 in (7.11). In terms of the new indices li, I2: 



II = 11' + qh - XN, X = iB.t[{n' + qli)/N] 
Ni=l2 + X 
n2 = ph - MI2 . 



(7.13) 



Also, the argument of the translation operator in (7.9) simplifies, since (using (7.10), (4.3) 
and (7.11)) 

n2 + —{ki- k[) = (ph - MI2) + —{Akh + t^k) 



K 



K 



/MAk 



+ P]h = —h . 



(7.14) 



After renaming some of the dummy indices, the Dirac bracket may be written in the form 

. 2 °o 00 00 00 
/(k',k) = ^ E E E E (-l)^^^^"^-''^^^+^^"^5(fc^-/c2-27riV2) 

JVi=— 00 N2 — — 00 ni = — 00 n2=— 00 



\N\-l 

X exp 



x5{ki — k[ — NiK, — 77-1 AA;) exp 
.fki + k[ 11 + 11' 



^2 + ^2 ^^^_MN^) 



-27ri 



2k 



n2 



{(j)^\f(^Ai + —niA2)\(p^+qm+\N) ■ (7.15) 

iV K 



Using the fact that the type II Wannier functions satisfy |(/)^_|_iv) = \4>ij.)-i the Dirac bracket 
(7.1) may finally be written in terms of a set of coefficients /nin2 the form 



00 00 



/I 2 °° 

/(k',k) = ^ E E E E (-1)'^''''+^'^'"'+™'^' <^(/i;2 - A;2 - 27riV2) 

JVi=— 00 N2——00 ni — — oo n2 = — 00 



xS{ki — k[ — NiK — ni A/c) exp 



X exp 

The coefficients Inin2 ^-^e given by 



n2 



•-n\n2 • 



(7.16) 



|Ar|-l 

^nin2 ~ ^ ^ exp 
/u=0 



27ri / 1 s 
— 2 9ni)n2 



.> • (7.17) 



8. Representations in terms of translation operators 

The expression (7.17) for the coefRcients defining the Dirac bracket, and the expression 
(6.8) for the Wannier function image of a translation operator acting upon a Bloch state, 
can both be expressed more elegantly by defining extensions of the magnetic translation 
group. 

First I will define a translation operator which acts upon the labels of the type II 
Wannier states. For integer values of Ai and A2 they are defined by 



i{ni,n2)\4>i^ 



exp 



27riM 
N 



(/X - ini)n2 



\4>^l-n^) 



(8.1) 



These operators were originally introduced in [7] . They have an algebra analogous to that 
of the magnetic translations: 



i(ni, ni)i(ni, 71-2) = exp 



27riM 
N 



{nin'2 - n2n[) 



i(ni + n[,n2 + n'^) 



(8.2) 



Using the definition (8.1), the coefficients Inm given by (7.17) which define the Dirac 
bracket (7.1) become 



liV|-l 



/n.n. = {-ir^^-^ J2 {<i>'M-Qnuqn2)f{^A, + -^A2)\<j>^) . (8.3) 

ix=0 



A further simplification can be introduced by using the notation |$) to represent the set 
of Wannier state vectors /U = 0, ..,|A| — 1}. The object |$) may be thought of 

as a state vector in an 'expanded' Hilbert space, with inner product 

|iV|-l 

($'!$)= (8-4) 

Equation (8.3) can now be reduced to a satisfyingly simple form by introducing a gener- 
alised magnetic translation operator in the expanded Hilbert space: 



t(R) = {-l)Pi^^-H{-qm, gn2)T(^Ai + ^ A2) , R = mAi + n2A2 



K, 



With this definition 



'■nin2 



($'|r(R)|$) 



(8.5) 
(8.6) 



Also, comparing with (5.9), it can be that the gauge transformation exp[ik.R*] results in 
a transformation of the vector of type II Wannier states of the form 



1$') = t(R*)|$) . 



(8.7) 



The operators T(R) again have a non-commuting algebra analogous to that of the magnetic 
translations: 

. (RAR') 
7ri7- 



r(R)r(R') = exp 



where 



7 



K 



(Ai AA2) 



T(R+R') 



(8.8) 



(8.9) 



is the dimensionless magnetic field parameter mentioned in the introduction (equation 
(1.5)). 

Prom (3.16), it is seen that the evaluation of the matrix elements of the Hamiltonian 
involves calculating the matrix elements {B{k')\T{r)\B(k)), where r = l3{i'iAi + V2A.2). 
The Dirac bracket (S'(k')|S(k)) was obtained in equations (7.16) and (7.17) in terms of a 
set of coefficients Inm-i • The calculation of section 6 shows that the operator T{y) acting on 
a Bloch state creates a new canonical Bloch state with k shifted to k + z/2ai. It is natural 
to expect that the Dirac bracket (i?'(k')|i?(k + z/2ai)) may be expressed in the form (7.16), 
with the coefficients /nin2 replaced by /$'$(ni, 77-2, ^'2) (note that /nin2 — -^'i>'4>(^i, 7^2, 0)). 
Noting that pn + MAk = 27r, and that ki is replaced by ki + 27^2 in (7.16), it is found 
that 

00 00 00 00 

(S'(k')|S(k+iy2ai))= I] I] J] ^ ^_-^.^N,n,+pN,m+MN,N, 

N1 — — 00 JV2 = — 00 ni=— CO n2=— 00 



x6{k2 -k'2- 27TN2)S{ki - k[ - Nm - mAk) exp 



k2 + k'2 



ipm - MNi) 



X exp 



n2 



/$'$(ni, 712,1^2) 



where 



/$'$(ni, n2, V2) = exp[i7r(p + 2tt / K)n2V2]Ini+Mu^, 



n2 



(8.10) 



(8.11) 



Now consider the evaluation of the matrix element (-B'(k')|T(r)|i?(k)). This may be writ- 
ten in the form (8.10), with the Wannier state |$) replaced by the state |^»") given by 
(6.8). Combining (6.8) and (8.11), the coefficients may be written in the form (8.10) with 
coefficients 



lA^I-l 



/$'^>"(ni, 77,2, J^i, 1^2) = exp[i7r(p+ 27r/K)n2J^2] ^ exp 



27riM/xz^i 
N 



K(3 



2txN(3 



x(0' |r(R + Mzy2Ai)T -r^viAi + '-V2A2 



K 



(8.12) 



This coefficient may be expressed in the form 



U'^"{ni,n2, 1^1,1^2) = ($'|r(R)'r(r)|$) 



(8.13) 



where 

T{r)=i{-iy2,iyi)f(^-^uiAi + piy2A2^ . (8.14) 

The operators T(r) commute with the T(R) operators: 

[f(r),t(R)] =0 (8.15) 

for all lattice vectors R and r//?. 

9. The generalised Peierls eflfective Hamiltonian 

9.1 A one dimensional effective Hamiltonian 

The motivation is to obtain an effective Hamiltonian, having a spectrum which is the 
same as a subset of the spectrum of the original Hamiltonian. The effective Hamiltonian 
is easier to analyse because the number of degrees of freedom has been reduced. The 
approach is analogous to that used in earlier work on the phase space lattice Hamiltonian 
[7,14]. The Hamiltonian will be reduced to a block diagonal form, and matrix elements of 
the Hamiltonian within one block will be compared with matrix elements of the effective 
Hamiltonian. If the basis states are in one to one correspondence and the matrix elements 
are equal, then the spectrum of the effective Hamiltonian is the same as that of the block 
of the full Hamiltonian. 

In the case under consideration, matrix elements of the Hamiltonian will be evaluated 
in the basis formed by a set of generalised Bloch states |i?'(k)). They are compared 
with matrix elements of an effective Hamiltonian Hproj in a suitable basis with elements 
\^{x,k2)), and the coefficients defining Hproj ^ire chosen such that the non-zero matrix 
elements of Hproj correspond with those of H, in that 

47r2 _ . _ 

(5(k')|iy|5(k)) = — (^(x',/c2)|iyproj|e(^,M)'^(^2-A;^) . (9.1) 

This equality holds when k2 = and where the states |^(a;, /C2)) are labelled by a contin- 
uous variable x — ki/n. 

It will be assumed that in the case where the dimensionless magnetic field /? takes the 
rational value p/q, there is a non-degenerate band. It will also be assumed that the gap 
separating this band from the rest of the spectrum does not close when is perturbed 
away from the rational value p/q. The effective Hamiltonian is constructed to reproduce 
that part of the full spectrum which evolves out this band when (3 is perturbed from the 
rational value. The type II Wannier functions |(/)^) for this band are determined, and used 
to generate a set of generalised Bloch states using (4.1). A projection operator P = f{H) 
is applied to these states, where the function f{E) is unity where E lies inside the band, 
and zero throughout the rest of the spectrum. The states resulting from applying this 
projection 

|S'(k))=P|S(k)) (9.2) 

are orthogonal to all eigenstates outside the band, and therefore represent the Hamiltonian 
in block diagonal form. The projection operator may be written in the form 



/CO 
dtfit) exp[im] 
-00 



(9.3) 



where f{t) is a Fourier transform of f{E). The stipulation that the spectrum has a gap on 
either side of the band ensures that f{E) can have arbitrarily many continuous derivatives, 
implying that this integral is nicely behaved. 

The projected Bloch generalised Bloch states are sufficiently numerous to form a 
complete but not overcomplete set for the band, and may be assumed to be complete 
provided the matrix clement {B' {k.')\B' {h)) is sufficiently small when k k'. This criterion 
can be tested and verified using the results of sections 7 and 8. Because the states are not 
orthonormal, a normalisation operator must also be calculated, such that 

{B'{k')\B'{k)) = ^(e>',/c2)|iVp.oj|e>,fc2))5(fc2 - k',) . (9.4) 

The subset of the spectrum of the full Hamiltonian which lies in the projected band can 
be determined exactly by solving the eigenvalue problem [-ffproj — ^'-^proj] IV') = 0) or 
alternatively by calculating the spectrum of the effective Hamiltonian operator 

^eff = iVp-oj'4rojiVp;i(' . (9.5) 

Consider the matrix elements of the Hamiltonian, expressed in the form (3.16), in 
the basis formed by the generalised Bloch states. The wavevectors k and k' can both 
be restricted to the first Brillouin zone, i.e. ki,k[ G [0, «:) and /c2,/s2 ^ [0;27r), because 
these states form a complete set. Alternatively, states in an extended Brillouin zone can 
be used, since they only differ by a phase factor from the states within the first Brillouin 
zone. States with ki differing by multiples of k are identical (apart from a phase factor). 
Similiarly, states with k2 differing by multiples of 27r are identical. When writing matrix 
elements of the Hamiltonian in a complete set of states, the summations over Ni and N2 
in (7.16) can therefore be dropped: 

(S(k')|i^|S(k)) = ^5(/c2-A;^) Yl S{ki-k[-mAk) 

ni=— 00 



X exp[ipk2ni] exp 



n2- 



-2m 



ki + k[ 
2k 



n2 



ttI 

nin2 



(9.6) 



In the case where (3 is rational, uk + mAk = for some choice of n and m. In particular, 
7 = Ak/ K is also a rational number, 7 = p' /q' , so that this relationship is satisfied when 
n is a multiple of q' . In this case, only q' distinct states are coupled, and the Hamiltonian 
is represented by a g' x q' matrix with parameters k2 G [0, 27r) and ki e [0,Ac/g'). In the 
general case there is no finite dimensional representation. 

Now compare the matrix elements (9.6) with matrix elements of an effective Hamil- 
tonian of the form 



^proj = i^proj(K) = ^if'(R) exp[iK.R] = 5^if'(R) T'(R) (9.7) 

R R 

where sum runs over all of the lattice vectors R = niAi + n2A2, and where the following 
relations hold: 

K= — (ai^i + a2^2) (9.8) 
[^1,^2] =27ri7 (9.9) 



(here the are reciprocal lattice vectors, satisfying a^-Aj = 27r5ij). The operators gi and 
§2 have a commutator which is analogous to the usual position and momentum operators. 
Eigenstates of §2 will be introduced, with eigenvalue x: ^21^(2^)) = x\^{x)). Evaluating 
the matrix elements of (9.7) in this basis leads to matrix elements which are very similar 
in structure to (9.6), if the coeflBcients H'CR) in (9.7) are identified with the coeflBcients 
H'^^^^ in (9.6). The correspondence becomes even closer if the states \i{x)) are 'gauge- 
transformed' as follows: 



\i{x,k2)) = exp 
The matrix elements are then 



. / pk2 



(9.10) 



00 00 



(e>',A;2)|4roj|e>,A;2))= E H'{R)exp[i{x + x')n2/2] 



ni = — 00 n2 = — 00 



X exp[ipniA;2](^(a; — x' — 2'K^n\) 



(9.11) 

are identical to the 



Identifying x = ki/n and 7 = Ak/n, these matrix elements of Hpro] 
elements (9.6) for all values of /c2. The spectrum of (9.7) is therefore identical to that of 
(9.6) when 7 = A/c/k and H'{R) 

Coefficients of the effective Hamiltonian 



TTl 

n\n2 ■ 



It remains to determine the coefficients -ff'(R) = H'^^^^ in (9.6). These are obtained 
using (7.16) and the notational devices introduced in section 8. The Hamiltonian is given 
by (3.16), and takes the form of a sum of magnetic translations of the form T(r), where r//3 
are lattice vectors. The action of the Hamiltonian (3.16) upon a Bloch state |-B(k)) may be 
represented in terms of the action of an image Hamiltonian upon the Wannier states that 
generate the Bloch states. The matrix elements of the Hamiltonian, (i?(k') |if|i?(k)), are 
of the form (8.10), with the coefficients /$'*(ni, 77-2, ^"2) replaced by coefficients H'^^^^ = 
if'(R) characterising the Hamiltonian. These are given by an expression analogous to 
(8.13): 

H'iK) = ($|t(R) n\^) . (9.12) 

The operators in (3.16) commute with the magnetic translations, and therefore commute 
with f(r) and 7'(R). Using (8.13) and (3.16), it is seen that the operator Ti, which is the 
image of the Hamiltonian in the Wannier function Hilbert space, is 



^ = E^kT(r(k)) 



(9.13) 



where r(k) = /3(n2Ai — niA2) corresponds to the reciprocal lattice vector k = niai+n2a2. 
The image H of Hamiltonian in the space of the Wannier states commutes with the image 
of the lattice translation operators: 



[7i:,t(R)]=0. 

A similar representation exists for the projection operator P = 

image in the form of an operator P acting upon the Wannier states. 
^ . . . 



(9.14) 

f{H)\ this has an 
Also, the operator 



proj 



VHV which is the image of the projected Hamiltonian i?proj acting on the Wannier 



functions may also be expressed in a form analogous to (9.13). The effective Hamiltonian 

can also be represented by an operator Tieff = V~^^'^'HV~^^'^ acting on the Wannier states. 

The formulae discussed above can be used to calculate the Fourier coefficients of 
the effective Hamiltonian using (9.12). Methods for calculating these coefficients as an 
expansion in (3—p/q are discussed in [14] for the case of the phase space lattice Hamiltonian, 
and these techniques may be adapted to the present problem. In order to establish the 
validity of the Peierls formula, it is necessary only to establish the limit of the coefficients 
if'(R) in the limit P ^ p/q. These coefficients are identified by noting that, upon setting 
(3 — p/q the Bloch states become eigenstates, so that {B(k.')\H\B{k)) — £{k)d{k — k'). 
The corresponding expression (9.11) for the matrix elements of the effective Hamiltonian 
reduces to 

oo oc 

{^{x',k2)\Heff\^ix,k2))= Yl Yl H'(R)6{x-x')eM^xn2]e^p[-ipn^k2] . (9.15) 

ni = — oo n2 = — oo 

In the limit /3 — > p/q, the coefficients H'{IV) of the effective Hamiltonian are therefore the 
Fourier coefficients of the dispersion relation. The effective Hamiltonian (9.7) is therefore 
of the 'Peierls substitution' form, (1.1). 
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